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ABSTRACT 


The applicability of the methods of system identifi- 
cation to the determination of ship model motion parameters 


is studied taking as an example a model of the Mariner 
class hull form. 


A mathematical expression describing the dynamical 
behaviour of an ocean vehicle is first presentede Then the 
values of the ship model motion parameters to be identi- 
fied are predicted using available techniques. Data 
representing the dynamical behaviour of tne hull form 
Mneey controlled conditions is obtained with model tests 
meme bowing tank. This data is then used together with 
Piet rst estimated values of the motion parameters to 
Semermire tne corrected values of those parameters. 


The sensitivity of convergence to several noise 
factors carecteristic of the process is also studied. 


The results indicate a fairly strong dependence on 
those noise factors and on the initial uncertainty that 
in general decide the final convergence. This fact seoms 
to indicate that the particular technique used and the 
Paericular input adopted are not in conjunction the most 
appropriate for this type of dynemical behavioure 
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Consider a system whose inputs are known and the cutruts 
observable. System identification, in general, has the 
Mijective of determining the structure of that system 
using the known inputs and measuring the outputs in order 
fie oie may be able to predict the behaviour of the 
System (outputs) to the various excitations expected from 
M@@emenvyironment (inputs). One such a system of interest 
an navel architecture is the ocean vehicle. The naval 
Meee 1vect wants to be able to predict the velocities, 
@emeerecions and other important motion resvonses that a 
mmo!) | exverimnent at sea in order to erfectively design 


mere orer response, structural integrity, and other 
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GTical aspects of operability 
The identification of a system of completely unknown 
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Pepuesyre iS, in most cases, 
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ossible, zis) | a 
being, There are, however, some types of systers for 
Beeeteeeentification is possible. Those consist of para- 
Mefereaentitication problems where the general structure 
of system is known but the particular values for 2 set of 


rs are to be determined. The ccsin vehicle was, 
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eevee tely, these characteristics, Using rigid hody 
Peeterecs, one is able to develop the so-called equations 


Semmocton wnich describe the dynamical behaviour of the 
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vehicle. In these equations there are some varameters 
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(ship motion parameters or hydrodynamic Coett ice eee 

are not known or at best whose value is known with 2 certain 
G@eegree of uncertainty, It will be necessary then, in order 

omsuliy predict the behaviour of ocean vehicles at sea, 


iomebbain the correct values of those parameters or, at 


( 


iMeast, oy sclected measurements, to decrease the dezree of 
Meer tainty. Unfortunately these measurerents are, more 
Sieeem Lian not, combined with noise from thea SWE er be 
bnemselves or from the environment et tne time the measure- 


ments take place or from both. These types of noise need 
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meee ae en inte account when identifying the pararetere, 
An examnle of the foreroing is the case of experiments in 
a towing tank to obdtain deta from models to anply to full 
Scale vehicles. Measurement noise exists whenever one 
Memeemony type of transducers that always chow some 
mereresiS and are characterized by a specific dymamical 
Benavyiour. Environment noise will be caused among others 
memeeerterence between the model and undesired waves, 
Siaevirration of test egquipnent., Of course these effects 
could ke included in the mathematicel model of the vehicle 
but it would make the analysis of its behaviour much more 


Sieorenilt. Thus it is seen that sometimes it’s preferable 
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to 2ssume aA simpler model for the hehnavicu sore spose 
and introduce the effects not accounted for in the mathe- 
matical model as noise (environmental noise) whose seneral 
nature 1s assumed known, 

Many of the available technigues of system identifi- 
cation were developed for a dynamical system whose 
behaviour may be expressed in a so-called state space 


representation[ 1]: 


% = Teas |) ; DAL las, oa aaa (ileg t) 
vy = h(x,u,t) | (1) 


Peteeemciec vector x expresses tne state of the system at 
meee (+> or) the vector v represents (eee > lou aol) 


a Mie Ra ai je > * we s + Jie any 
memeeac time = and the vector y ns 14017 25 1 ee. 


me as the state of the system at time ice Gear olen dr to | 


and x is the derivative of x with respect to time. f and 
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meere vectors that describe the structure of the system. 
TRE main assummtion to be made at this stage about the 
meemensystem in crder that meaningful results can he 


Merete co is that there will be wnicue values of =«(%) and 


Meemece = Siven u(t) and x(t.) = x, for all t such that 
ee, 
Sis 
Mmeepurpose of any technique of system identi tieation 


femeeinen be, by opernting on u(t) and v(t), the deter- 


fiweeesmon Of f and h. As said tefore it will b*, in 


1a tun 











eneral, not possible to obtain the solution vo sie sure. 
Therefore, the assumption of a known structure 
except for a set of unknown parameters is made and the 
system is expressed by tne following: 


oe ane Gn Gene omer, x(t,) =x 
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Meee Le VeCtor p represents the set of parametcrs te 
be determined. 

it: order that the identification of those parameters 
Eememescome value, it is further assumed that p is time 
Beeeeriant in the time interval of interest t,< tt, or 


ewes, Ce tS «6Utiine devendence is known mesning tnat 
~ | 


once some n(0) , Ue Sei , is determined, De) widd 
Pemeievn in the interval [tos eae 

iemust be emphasized at this stage that the fore— 
going does not assume the unicueness of that representa- 


tion, Only the uniaveness of the input/output relation- 


assumed, That is, there will be for a given 
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Peuee tb) and 
= 0 

y(t) = blul(t 
mmo be possible to snow the existence of a repre- 


pemereeron x", f*, h* and p* such that 
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kt = f£'(z',u.0'4¢) Ta G coe 2 LS) 


y = h'(x'u,p',t) Be, 
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The problem is then redefined as the determination 
Oi a set of parameters p 02 Some Stabe wepeecon | teem mes 
the system defined by y(t) = blu(t,,t), +t). 

Comins back Wer theo ewer an ocean Vente le, a tweed 
Beeshovn that the resulting equations ef motion may be: 
resardcd as a state representetion of the vehicle 
G@yneamical behaviour and that the assumptions made so far 
are met by that mathematical model, 

The foregoing assumes, arart from the unknown set 
Seema cameters p, no other tynes of uncertainty, But, as 


said before, it is sometimes advantageous to select a 


pla 


ieeeer meatheratical el-for @ Siven. Sys bem and teo- 
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evce t! alled system or environmentel neise since it 


1@ ¢o0-c 
mescwies irom interactions between the system and the 
environment, Also consideration must be given to 
measurement noise. A model that takes into account those 


Bateters sill be of the fellowing form [eu : 


xd = ieee, Ds v) x(t.) = x (138) 
y= h(x.u,d-¥,¢) (1,9 


Peeors in the structure of the system will come through 


Meme ironnental noise. 
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In addition (case of time varying systems) the parameter 
Pee Lon Dp Nay have the 1oOLLOWwIng sot ruietmmee. 


OS Ge a) Dit eo leo.) 


where £ 1s known and w’ is noise, 

In order to simplify the above formuléetion the 
following assumptions are generally made: 

lt, The mathematical model 1s time invariant. 

2. soyetem structural uncertainty will appear as 
an apparent additive noise process, 

3, Measurenent noise may also be considered an 
additive noise process, 


Wen ate pale oom 
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ry 
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me Ouvput roasy 
iivesstvete of the system and are structurally indevencsnt 
eieetre Input, 

5. Parameters are coefficients in the medel 
[ewe w7re that are to be identified. 


Tne system dynamical behaviour may he expressed 


x = f(x,u,p) + w Gib) 
2 ae (1,12) 
mei kt Vv ee) 
Wiere is the measurement matrix, often a diagonal 


Marx. 
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The purpose of the identification will be the deter- 
mim2bion of Dp. 

In Chapter 2 the ecuations of motion Tor am ecean 
vehicle will be developed and it will be shown that indeed 
they constitute a state space representation of the ocean 
vehicie’s dynamical behaviour. It finial ig be shown also that 
they satisfy the foregoing assumptions. 

In Chanter 3 the methods of predicting the ship 
model motion parameters will be outlined and values of 
meee parameters for the mariner class hull form will be 
Mieeeented (Anpendices 1, 2, 3). 
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Will be described and its several components isolated, 
In Chanter 5 the experimental procedur?2 followed 
a7order to obtain the data will be described, 
iMomonapter 6 the different methods of system 
Meemeerication will be discussed and those which seem 
more convenient will be used to identify the motion 
Bemeeme ters of the mariner class hull form. 


eo? the- results bce imodawe Opa ter soe ie 
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In Chante 
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@emeerecussec and in Chapter @ the consequent conclusions 


mieeeeconmendations will be pointed out. 
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eoUATION OF MOTION 


In order to develop the equations of motion for an 
Mees) vehicle, it 18S convenient first 6 @e7 mie wee 
rectanfular right handed coordinate systems, one fixed 
Meeene surface of undisturbed water with the positive 


axis pointing downward and the serond moving 


fe ~«C“o 


merece tne vehicle, the origin of such ecordinate system 


° ° a * pe a De = \ 
being arbitrary (but the coordinstesc (fail er 2a) 


meme, Of eravity of the vehicle in this system boing 


known) pointing the positive x, = x axis to the bow of 
m@eemvyenicle and the ee Pointing downward in the 


Sear labrium position (see Figure 2,1) 
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The orientation of the moving 2xis system witaeses 
to the fixed axis system will be expressed by three angles 
(the “modified Evler angles") defined as follows, 


assuming that initially both axis systems coincide [3]: 


ae 


a angle of rotation of the movingwetteee ee 
about the z = z, axis (yaw angle) 
@ - rotation of the moving coordinate system already 


yawed about its v axis (pitch angle) 
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otation of the moving coordinate system 


feieesady yawed and vitched) about its x axis (roll angie). 
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Defining the rotation about the body axis 2s 
Peepectivel: §,(about the body x; axis), Ss(about the 
body Xo avis) end Sg(about the body x3 axial tae Yeleaia. 


between those rotations end the Enler angles will te (3); 


eostsint+ 3cos C25, 


Mae motion of a rigid body in six desrees of Treacom 
Peeemeresnect to an axis system fixed in enree, when the 


Seetie is at the center of gravity and the moving axes 
Me@eere Principal axes of inertia, may be represented by 
Deron S equitions, nanely: 


(ru) 


fee 1S 2 diagonal mizrix whose elements are eguval to 
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the mass of the body, Up LS the translational are oe aa 


mector of the center of fravity, of components Higa =a. 


Moving axes u_,; Vq and War Eis the total vector force 


G 
Eeing on the body of components along theraxycsu* 7. 1s a 


@, 1 is a diagonal matrix whose elements are the inertia, 
Beis the rotational velocity vector ehout those axes 
z 


Beecomponents p oa uate! ae and M is the (totale@ier. 


Caer 


meer On the body of components K, M, and N, Negleeting 
Meewetiects due to rotation of the carth in ocean 

vehicle dynamics Abkowitz [41] and Tufts [5] develop 
Memations (2.2) respectively with fixed and verving mass 
Mmeemposition cf center of gravity for the case of an 
ee ary oricin, the firet vsirng vector enelysis and {he 


second Lacrang?'s equations, Therefore only the resulting 


equations will be presented here and for a body of fixed 
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if 


es: ae 
mi utqw-rv-x,(q tr )+¥.(pa-r)+z.(pr+a) | 
- AZ Z fi 


6 ? 2 r e os 
ie = ml v+ru-pw-y)a(r +) yz, (dr-p) ex, (Cee (Ze? 
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Z ; ; 
Z = m{ftpv-au-2 (p+ )4y(rp-4) ty, (rat) | 
; : Ce Coes a! 


° ° e ; 2 ? 
= TP tie (1-P2) + ee Cane) < To 64 aaa ) — 
(Tt yy) Germ, ee peten) 1: my (w-Guipv) oe 
Meet ag +i ({r-qp) + 1I.,‘prar) + I Cay + - 
ty ye ye es 1 an 
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Beoere unsubscrioted velocities 


ep) 


ind eecelerationse, 2 7eo 
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imeean arbitrary origin and x 


a) 


Ye and Z. are che 


G’ & 


meerdinates of the center of gravity. 

Note that F and M will in general ve functions of 
position, velocity and acceleration of the yehicle as well 
@eeproperties of fluid, geometry of the body and of the 
inouts. For eae body of given geometry and mass distribution 


ana in 2 siven fluid then 
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merce the purpose of this investigetion mes tne 
etudy of the aprlicability of system identification 
Meemmaques to ocean vehicle dynamics and not the 
identification of parancters of a particular vehicle 
(re 
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ameters imbedded in F and M™) it was decided at this 
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peeety the formule tion Dy Peon eae = 
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Model to one or 2t mest two dogroees of freedom. It wall 
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P eens vderod heave (traneiat en 2 lems 6a eee is) 
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BEG # : : ° ° 
10207 ohout the ~ ~yis) with genre vorseand 


speed. To simplify the formulation even more eee etl 7 
eal 
were mane zero, Merci ore, the equations of motion 


Meee this restrained mode (y,=2 =QSurvaveaq=G=Pr=r=x =y =e) 
will be 
te =. HY 
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Fieure 2.2 
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Since Y, 9 and Vv are zero, the 2dove relate 


medyce to 


Os9 (2.8) 
Bigs = W cos D a 


Ignoring for the moment the presence of w and bd as 
arguments of Z and K the complete set of equations in 


state space form would be: 


ee L(y WZ eDa PrP inout) 
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D = a K(% 4.2 De PoP input) 
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iemmerder Oo ne able to vrocesd end aoply the existing 
identification technicues to ocean vehicles (specifically 
Pieeree Shics}, it is necessary to isolate w and that 


the dependence of Z@ and K on w and vo is linear. Therefore, 


iScume thevexi1= ter eeser 


4. - 


a 
o) 
oF 
e) 
99 


the next stern will 
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Mavilor’s expansion zbout 
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T 3! Bwbw Di ewoe | 200" ou are ened 
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and in the same way for K. 


It 1s known from potential theory [5] that for a 
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Tae ae 
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bedy moving in an luid away Tron the surfa 
Miemo are no dependence of force or moments on products 
Seececelerations with veloc ities or accelerations. 
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Sieeeace vehicites without any other justification, and 
Bieeeees also all other terms where acceleration is 
multiplying any other variable. 
Mocr the vehicle in equilibrium also Z(o) =0, 
Purther simvlificetion can he achieved if we 


Scroider eich element of tre expansion, which will be 
considered in the nert chapter. It will ve enough to 
femeeo. the moment thet due to port and eisrboard symmetry 
Mem meoonstrically and inertial) the vertical motion 

of the vehicle does not produce any rolling moment, at 
Messe to first order, and therefore those terms car be 
Meee tron the equations furth 


Pee yy, and to first order only, rolling motions will 


Meme OOcCe any vertical force or, in otrer words: 


Bee OM 8K eK | esi) 
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Also the square terms 6*2 and 8&?K will be nonexistent, 
ow~ 60" 


instead, consideration must be given to Re iepae need 
moctical effects which are of the form cyw|v| and cP |p| ° 
Although these terms are not analytic, they will be 
miewuded so that the model represents a closer approxi- 
mation to the vehicle physical behaviour. 


The form of the equations will then be; 


Maoh [Awe pw dwltZ, tote et 1 nyo Gy) 
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foephe next chavter each term in 
memes Of first order, quadratic or at most cCubie will 
be discussed in relation to its nature, order of magnitude 
Peewee nods of predicting its Galles. 

ieecummary, from Newton’s equations and for a case 


Of restricted motions, assuming the existence of a Taylor 


Series for the force and moment and Linear dependence of 
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mec lON PARAMETERS 

For convenitonee of presentation, the escvce uae. 
mee motion naremeters will be subdiviced into threacmer, 
S-oups, namely: 

m) Motion parameters resultant of @ifterentua me 


Seeeeorce or moment on the vody with resvect to position 


St 


(7,4) only (rydrostatic coefficients 


me) linesr terms resultant from differentiation von 
Memee or moment on the body with respect to acceleration 


These coerricients are dependent only on the geometry 
Beomease Aistritution of the hody, for example, share 
Br@eretetive vorition of the metecenter, center of 
Merete, and centor of gravity, for the restoring moments, 


Memeo oetricients sre calcvleated based on the Vaaec 
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meee or the hull form, 


Meme te Shim in tre upright position, the Myevent 
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eelected value of Za) Then at each Waterline aimee 

of roll is simulated in the body plan, being the point 

eo cotation the assumed position of the center of sramaays 

and the resultant buoyant force 8B (zn ® ) and the fign tan 

ar”) GZ(25,9) are calculated (for details see, for examole, 
Beeanciples of Naval Architecture" edited by J. P. Comstock, 
Miepver 2). The righting moment ( -K(2, 5) ) is then 

given by the product of Bae?) , CZ( ZS The vertical force, 
 meSsitive downwards, is the difference between the weignt 

meene body and the buoyant force, or in other words, 


Mensiderings the comoonent of the force along the body 


~ “ ©@ 


22049) = (W ~ R( 2546) 206 p 
K(25,6) = -B(zysP) + G2(%99) 


wnere W is the weight of the vehicle, 


feoeorasc, ie fave: 
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ime values of Z2(: 
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0) and K( 254) are usually 


presented tabulated or in graphical form (see Apperdix 1). 
feemeeder to get the hydrostetic coefficients, Aerivatives 


Memes), it is necessary to obtain a mathematical 
representation of those ee for example, a pelyromial 
Mmeemere bY the method of least scu2res. Cne such se nere 
was develoned by Kerwin [4] end implemented by Parissis [7] 
Beene Lerendre Polynomials to represent ship forms whsch 


Peeedually applicable to (3.1). 
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CHAPTER 6 


MOENTIFICATLON OF SHIP MOTION MODEL PARAMETERS 

Most of the available techniques of parameter identi- 
fication have been developed for linear dynamic systems 
[1]. Some of those techniques, assuming the linearized 
equations of motion to be a valid description of the 
behaviour of an ocean vehicle, have already been applied 
in the determination of ship motion parameters in heave 
and pitch namely, the frequency response [27] and the 
step response methods [28]. 

In the case of a non-linear model, two approaches 
are open, both yielding approximate results. 

The first, making no assumptions about the character- 
istic of the uncertainty, assumes a mathematical model 
for the system and then comparing the output of both 
model and systems to the same input, searches for the 
structure of the model that minimizes a function of the 
error between the two outputs. An example of the fore- 


going is the model reference technique, For the system: 
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where u vector input to both system and model 
Bite = vector output of the system 
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f, - vector structure of the system 
_ - vector structure of the model 
i - scalar function to be minimized. 
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where the values of Dp. are unknown and the values of 1B 
Will be determined by minimizing V, which is often of 
mie form 
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ee Po <t and may not apply outside tnat interval. Also 
it may not be true for non-linear systems that 


‘La se = Go) (6,4) 
t—" 00 m S 


as one might expect from linear systems. The optimum 
value of p in Equation (6.3) could be found by a maximum 
m 


gradient seeking algorithm (first order gradient method) 


or by a second order gradient method [| 29] 
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The second approach is to use linear estimation 
techniques, such as a Kalman filter to the linearized 
form of a modified state representation of the system. 
Considering a linear system driven by noise, Bryson 
and Ho [29, Chp. 11 and 12] present the following expressions 


for the estimate of the system: 


i) Linear Multestace Systems 


Bea 7 Fee * Sas 
x - state vector ELx | = x 

Blx -¥ ] Cx -x ]" = M, (6.5) 
w - environment noise ELw] = w 


a Jk 
ee ae w.-W, = O2b 
Cw, By] Cuj-w,T? = 0585, 
elw.-w.] [x,-k ] = 0 
Le —O =O 
4 G, [eens i GLon Matrices 
eee? = HW. X. + iV, 
ae 1 i TH 
fee measurement vector 
H; - measurement matrix eietey) 
V, - measurement noise ELy, | = 0 
vik 
™ Vv ; Pas! ae 
ELy.vs] = RSs 
-- gh 
BLX =x | vy =0 
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Assuming the sequences i of the different vectors 
to represent Gauss-HMarkov random sequences, they present 


Mme solution for the estimate of the state x —=—— * )aca 


X, =X, + Ki(z2, - Hyx;) 
— A 
ee hy Xe tT Og WG 
= Mv a Pa 
Ko = BHD Ry (6.7) 
. 21] SA? a ae 
P.=(M, +H, R; Hy] = 
a us ard a -l a 
My; - HyH, CUM Hi" + Ri) ~ Hews 
i _ F.P.F.) + 4.9.6, 
“$41 ~itisi 4i=4i4i 


where P, PouNererror cOVariance Matrix of Xs after the 


merenrement wnile M. is the error covariance matrix 


BG 
before the measurement. Similarly sa is the estimate of 
the system after the measurement while Xe is the estimate 
of the system before measurement, 
2) Linear Continuous Systems with Continuous 
Measurements 
The mathematical expression for this type of 
system is obtained by letting the time between stages go 
mem2zeroO, Therefore, for the system 
x = Fx + Gw 


ate), 


Sk PV 





with BS am. FCG ee = 0 
2 j } — 

G = lim Ga 

BO) aan 


Under the assumption that x represents a Gauss-Narkov 


randon process, the expression for the Kalman filter 


will be 
$= phe cw + K(z - HX) 
K=PH Ro 
P=FP+PpFl+cgo’-pH Ri uP 
ie (t,) == P 
inemefor the system as specified in (1.h)) to (i513) 
x = £ (xju,p') + w! (1,11) 
p = 0 GLI) 
Z = wx az AL es) 


the above technique is to be used, it is necessary by 


mmest to form the extended state vector 


a «| X 
Cemloe 
PD 
The new state equations are as follows: 
m= f(x,u)+w Elx(t,)] = x, ) 
Consleis) 
Z=Hx+V 
Miere in addition to 6.10 
a) |W" 
cone) 


0 
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JH 


fs | 
| (6.12) 


= |r] 


Now the state equations are linearized in order to 


lO 


i= 
[ao 


mee Equation (6.9) ors 


me f(x,u) + w + K(z-H x);x(t,) = x 


— O 
K = P HR” (6.13) 
ia r t-1 - 
Meomor P+ P $f +Q-P HR HP; P(t.) = i. 
ie ae 


The next simplification will be, considering that in 
practice most measurements are discrete in time, to 
approximate the above differential equations by a set of 
difference equations and use the formulation for a linear 


multistage processes or,in other words, for the system 


below : 
ey Xi + A-2h%; uy) + ArH; (6.14) 
ee 2 xX; + 


iime resulting filter will be, assuming for convenience 
that 


ELA *w; | = ou 
A ~— peak 
Key 7H, + acLOsi uy) + oKi fz; - Hx) (6.15) 


eee (R.A) 
K, = PH (R 
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Zh a -1 
= fit, = eo}. uw. ; : 
PB, = My - AMAH (HM, Hat 8-A)" HOM, 
(G5us) 
=i =i 


where R and Q are constant matrices, If the value of A 
is sufficiently small, this formulation will represent 
closely Equations (6.13) and will be much easier compu- 
facvionaly. 

In the above formulations, the partial derivatives 
may be evaluated in a nominal path or if greater accuracy 
is desired they are evaluated at x = x 

Bryson and Ho [29] mention several authors who have 
achieved success using Equations (6.13)and (6.15) for 
non-linear systems but at the same time make the warning 
that convergence may not be obtained if the initial guess 
is poor or if the disturbances are so large that the 
linearization is inadequate to describe the system, 

Moom the foregoing Equations (6.2) and (6.3), 
Pouations (6.13) and Peneene (6.15) it was decided to 
use Equations (6,15) due mainly to computational 
convenience. Due also to the poorness of the available 
data, as mentioned in Chapter 5, it was decided to atandon 


the identification in heave and coupled heave and roll 
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and concentrate exclusively on the identification of roll 
motion parameters whose data at a frequency of excitation 
of .66 cycles/second (roll natural undamped frequency) 
seemed reasonably good. 

As mentioned in Chapter 4, Equation (4.17) the state 
equations for roll motion, ignoring for simplicity the 


dependency of the system on the derivative of the input 


el © § 
p= K, p+ [ Kg - W.3BG ] (p-4)+ 
Jae ea aS. 
W. + DNS + [1/31K > + 
eee + (1/21, | Pipe L1/3tK * (417) 
p XxX <K°+T -Ke+I 
10) pe ja eo 


1/3!W.BG | (ely + -[1/3!W.BG | o° 
-Ki+l z =e oo 2 
Dp XxX Dp “xx 
P= p 
Non-dimensionalizing the above equations with the 


frequency of excitation the following is obtained: 

F = im 0 2 Al K ee cee 
Medea ee. Poot Sa ie 2 ee 

2 K+ wo | | w Re ae Sone 


| (ome) 
[p-f]+ | Wee 1 4 * 1/2 1K tpl}? “| P| 


Seeieatiints altel 


ee < 
olen ene W 


= eee = 


ARIEL 





-~ 


1: Bee ‘ 


(e]-[4 


Based 


2 


S77: 


IBA NeRG  Lalies =j 17s Peele 
-K*+I s( Sy For 14 
(OAS), 


Ae: 





pe xx XX 


on results of Appendices 1, 2 and 3, the 


following non-dimensional parameters will be,for 


w= 4,145 sec. 71; 


Ky s 
Katt U) 
K 1 


“Kitty uw)? 


_ (rr 
= e 2 
Bt tex GQ) 


iy 2iKk 
/ Pp Pp 


“Karly 





x 


31K ~ 


“Kotlyey 


1 
uw” 


miso and to 


between +1, 


SP 


ElO-. Elo 


~@ 0042 


= = a ee ( 30) 
= -,989 = x (4) 
= 2,744 = LOMedc> Weoley 
= les: = Be) 
es el a en 0) el 72), 
have the values of all the state variables 
it was chosen that 
= oy), 
= 1 x(1) C6 18 
ely (2) 
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ee 2), 
ious ) een 
The complete set will then have the form ;: 

ely(1) = .O1x(3).1x(1) + [x(4)-10x(5)] [.1x(2)-.4u(1)] + 
dO (5) to G2) at ; 
+ x(6).1x(1).1|x(1)| - [10x(7)-.167 10x(5)| 
(.1x(2)-.1u(1)] °= .167 10x(5) [s1x(2y ° 

miv (2) = .1x(1) 


or, rearranging 


y(1) = .01 x(3) x(1) + [x(4) - 10 x(5)][x(2) - uli) 
Pe ee 5) x (2) ae (6.19) 
+ «1 x(6) x(1) ABS x(1)+[e2 x(7) - .0167 x(5)} 
- [x(e2) “ u(1)]° - 0167 x(5) x(2)° 
cate) = x 2) 


Before using the data, this model was tested by si- 
milating noisy measurements —= only x(1) and x(2}) were 
measurable in the simulation -—- and the control Steal) 
with a format similar to the ons obtained in the data 
Collected. 

Also simplified forms of the state equations were 
Wesved. Their expressicn were : 

Non-lincar equations of motion to Identity linear terms 
mal) = «OL x(3)- x(1) > [x(4) ~- 10 x(5)|[x(2) - w(1)] 4 
elcome 05) oe 2) ee SOLES eee me | ieee ume 
7 2 (6.20) 
+ [.03123 ~ .0167 x(5) |[x(2) ~ u(1)]° + 


Wore x(5) <(2)7 








a [a 


C2 =x GAD) (cy2o) 


/menear Fovations of Motion 


me) =" 01 (3) ex Ger x5) =10 (Sie e 
blab) @| aie oxi S)) oe CORE) 
mec) = x(1) 


Pauarlvonses 6.19) 4006-20), and (6.21) were epeciwed 
varying the following parameters: 

im Initial State Estimate —- The “trie values” of 
the parameters (x(3) through x(7))were set in the simulation 
Mmeene values given by Equation (6.17). The initial 
Setimate was set either zero or Equation (6.17). 

2. Initial State Estimate Error Covariance P, - This 
is a diagonal matrix of equal elements which were set to 
Mec. and 5. 

3. Order of Magnitude of the Elements of R Relative 
to the Elements of Q@ - Both thes matrices are diagonal and 
of equal elements. The ratio of the elements of R to the 
elements of @ was set atl, 2, 5 and 10. 

4, Noise Level in Environmental Noise - w and in 
Measurement Noise - y ; Jo. a8 Osan sal 

Then the collected data was used in the identification 
procedure introducing small variations in the parameters 


mentioned in 1, 2 and 3 above around .a reasonable set 
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obtained in the simulation to obtain the best possible 
convergence and also to test the effect of initial 
conditions - or in other words, the necessary accuracy 
in the theoretical prediction of the values of the para- 
meters, to obtain convergence. Some of the results from 
the identification are presented in Appendix 5 while the 
discussion of the results is presented in Chapter 7. 
Summarizing in this chapter some of the available 
techniques for parameter identification were discussed, 
Attention was focused in those that permit application 
to non-linear dynamic systems. From those, one was chosen 
that having perhaps less accuracy in the results,is easy 
to implement. The equations of motion for roll went 
then through the necessary modifications to be in agreement 
feepiethne format of such formulation. The structure of 
the model and the formulation were tested by simulating 
the excitation and the measurements, and varying selected 
parameters to obtain a better insight of what to expect 


Mmeom tne identification, 
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CHAPTER 7 


DISCUSSION OF RESULTS 
1. Simulation 

The first equation to be tested was Equation 
6.19. The noise level was set to zero, the ratio of ele- 
ments of R to elements of Q, from now on expressed formally 
as r/q by similarity with the scalar case, equal tol, 
the initial state estimate error covariance matrix Pp at 
[2] and the initial state estimate x, at zero (the "true 
value" was set in all simulations at the values expressed 
in equation (6-17) ). The values of the astimate blew up 
and after 200 measurements the computer indicated overflow. 
This type of behaviour with no noise decided the use of 
Equations (6.20) or (6.21). These two models were both 
tested at no noisee Equations (6.20) wore first tested 
with x(0) = Eq.(6.17), r/q = 1 and P,= [5] The error on 
the estimate (XHAT ~ XX) where XHAT is the estimate and 
XX is the "true value” of the state,remained at zero during 
the whole experiment. However P(3,3) (variance of x(3) - 
linear damping coefficient) kept increasing with small 
oscillations. The other elements of P decreased quickly 
to less than .l1 . The same behaviour was observed for 


Equations (6.21), showing that the influence of the 
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non-linear terms is small, at least for the no noise 
condition. 
The last experiment in the no noise series was con- 


ducted with Equations (6.21) P 


P= antes r/q = 1 and 


(XHAT) 5= [Ole 

The error (XHAT - XX) for all variables converged 
rapidly to zero - in less tnan 100 measurements with the 
exceptgon of (XHAT(3) - XX(3)) whose value remained 
constant throughout. Again the term P(3,3) aeoaetels 
increased; while all the other variances came to values, 

after 500 measurements, less than .060 (initial value was 
2.0) the final value of P(3,3) was 2.202. The foregoing 
Suggests that when no noise is present the system converges 
fee nde ntiy or inivlal condiGions, with vie excem.uren 

of Equations (6.19) and with the exception of the linear 
damping coefficient x(3) whose identification seems very 
difficult to obtain. 

The rest of the experiments were performed with noise 
tevel equal to 1 which is believed to be quite a hign 
level due to the scattering observed in the megasurement 
vector Ze 

The first simulations weve taken with r/q= l 
XHAT>9= [OJand P.= [3] and = [5]. Both simulations blew 


20 
up at about the 140th measurement. At that time 9(3,3), 
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P(4,4) and P(5,5) were quite large, but different from 
one case to the other. The main point to be observed here 
is the conjugate effect of noise and not accurate initial 
estimate. 

The next experiment put together Equations (6.20) and 
(6.21). For both systems r/q= 1, P,= [0] and (XHAT). = | 
(XX) 4° Since initially the error covariance matrix was 
set to zero, it increased and then oscillated about 
Values equivalent to the ones found in the no noise 
condition (less tnan .950), with the exception of P(3,3) 
which steadily increased to a final value at the 500th 
measurement of .249 in both cases, with tendency to 
increase. The errors on the variables oscillates around 
zero with the exception of (XHAT(3) - XX(3)) which remai- 
ned constant at about .055. Again it seems that the 
identification of x(3) will be dificult to obtain. 

Next an experiment was performed with Eq. (6.20) 
with r/q=1, (XHAT), = (XX), but P - [5] which diver- 
Boa again confirming the influence of tne initial con-- 
ditions when combined with noise. 
The next set of experiments were performed with 
the noise level at 1, as before, but r/q = 2, 5 and 
10. 


Comparing the results from Equations (6.21) for 
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r/q 1 and r/q 2, (XHAT), = [0] and P, =[2, in the later 
case the estimate does not blow up although the error 
on the estimate and the variance both increase for x(3) 
and x(4), while for x(5) the error on the estimate and 
the variance decreases, which demonstrates some impro- 
vement. The cases r/q = 5 rfq =10 represent no advan- 
tage over r/fa= 2, presenting the same trends. 

Similar simulation was performed with Equations 
(6.20). Abain an improvement is noted when r/q goes 
up to two but no further improverent is obtained when 
r/q= 5 or 10. On the contrary, the‘decreasing in the< 
error seems to take place more slowly in the later 
cases In the case r/q — 2 the advantage of considering 
the nga REE terms when identifying the linear para- 
meters {3 noted because P(3,3), which increased in the 
case of Equations (6.21), is seen to oscilate between 
approximately 2.02 and 2.16, and P(4,4), after increa- 
sing for the first couple of measurement, markedly 
decreases. The Senco results were then compared 
with the case of P= [5], with clear advantage for tho 
POrme Yr » 

The foregoing seems to indicate trat the values 
r/q= 2 and Pj= [2] are the best suited among all 
miner combinations testéde At this stage small varia-= 
tions could be tried to further improve the convergence. 

But since the caracteristics of the data are not 


expected to follow very closely the caracteristics of 
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the simulation, it seems a better idea to try to iden- 
tify the linear parameters with the data obtained from 
the experiments and then vary the values of R, Q and P, 
in order to obtain a better convergence. But before 
that one last simulation was tried for the same values 
of Q, Rand P,, but for values of (XHAT), closer to the 
true values : (XHAT) = Keo oe, Sass ects | eivaue 
1/2 of the true values ). The number of measurements 
was also increased to 1000. It was noted that : 
i) P(S,5) decreased rapidly to a value of 4 
(after 10 measurements ) and then decreased oscillating 
until it reached a value of .057. The error on the esti- 
mate, after 1000 measurements, was oscillating between 
Bee 6 
ti) P(4,4) decreased less rapidly to a value of 
0566 ( after 30 measurements ) and then jumped to a 
valuo of .06 where it kept oscillating for about 350 
measurements decreasing finally to a value of .025 where 
it stayed oscillating until the end. The error on the 
estimate oscillated between +1 for 400 measurements 
and then decreased to zero oscillating about that 
point. 
141) P(G,3) increased to a value of 2ee2 after 
400 measurements, oscillated for a while and finally 
decreased to reach at the end a value of 1.29, with 


tendency for further decrease. Ths error on the estimate 
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increased to a value of 3.5 after 420 measurements, 
oscillated about that point for a while and finally 
Gecreased, osciileting to a value off sieat thesenda. 
Note that the turnover point seems to be the 
400th measurement and that a set of initial conditions 
on the estimate closer to the true value really impro- 
ves convergence. This fact confirms the believe that, 
Sentrary to the linear case, the initial conditions on 
both the estimate and its covariance play an important 
m7 the success or failure of the identification. 
2. Identification Using Experimental Data 
Following the conclusions of last section 

Pern was made with Q=[.1] , R=|.2] and ee [2%] ue 
Mae set of initial conditions was as prescribed by 
Equations (6.17). Simulteneously, runs with the values 
of Q, R and P| as indicated below were also made. The 
set of initial conditions, however, was not changed at 
this stage. It seemed preferible, since the above set 
is believed to be, among the available predictions, the 
elosest to the true value, to find first the best set 
of Q, R, and P, changing (XHAT), afterwards. 


Values of Q, R and P, : 


Q= [2] | 
R= [21]. 62], [4 
Po= [i .] ; [2.] : [5.] 


in the following combinations of R and P, : 








-~fl = 


MNGoS) sy) 

(fr ].f2J),¢ fry [s] 

( Fe ].feJ). ¢ bel].fed). ¢ fe] [5] 3 
( [-3],[1.] ) 


For the group ( esl iees ) the results obtained 
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with the first measurements were more or less as 
expected. P(5,5) decreased rapidly, P(4,4) decreased 
less rapidly and P(3,3) increased steadily. Then P(4,4) 
came close to zero and became negative for a while, then 
crossed zero again and went up to 2.48, decreased to 
ee/O and increased again to the final value of .42]. 
Beco, at che end P(5,5) = «#059 and P(3,3) = a475 .« 
Somes other runs a similar behaviour was noted. For 
cae. 2 | and ie [ce » P(4,4) crossed zero and kept 
decreasing up to a final value of -1.589. P(5,5) 
decreased to reach the value of .064 ( at the sane 
point at which P(4,4) became negative ) and then 
maeroased to reach a value of .123 at the erds P(3,3) 
increased steadily to reach 1.462 at the end. For the 
case R=|.2], Po=[5-] it was even worse since at 
peomend P(3,3) = =-4368225 and XHAT(S) = 252, /0.jeror 
R—|-l|, P,=|5-| 5 P(3,5) increased steadily, P(4,4) 
G@ecreased to a value of .175 and then increased again 
momreach at the end 425 with XHAT =-14.513 ( ? ). 


P(5,5) decreased to a value of .028 and kept oscilla = 





ting about that value, but with XHAT(5) always increa- 
sing. For R=-|[.1], P= [2] the behaviour was similar 
but not so extreme due to the smaller value of sae 2 
For R= [|.05] , P,=[2.-] , P(3,3) always increased 
while P(4,4) and P(5,5) decreased first to low values 
( .60 and .20 ) and increased afterwards, never becoming 
negative. For the lest run of this group ( R=[.3], 
Eee | 1. | ), the type of behaviour was similar to 
( ge -2 |; P=(1.] ) but- much more extreme —— the 
last recorded value of P(4,4) was -226.1 at the 350th 
measurement because some cycles later the computer 
stopped by overflow. 

The fact that P(4,4) became negative posed the 
meestilon of the correctness of the computerized schemes. 
However the program after thourough checking have been 

working successfully with linear systems. Besides no 
such strange behaviour was found during the simulation. 

Therefore the program was checked again in detail 
but nothing could be found able to explain the foreen* 
GOiINng. 

It was sugested then that as P(4,4) approaches 
aero, the computer agiatd) Hames some errors thus 
making it negative. The process would be selfsustaining 
afterwards. Assuming that to be the cause, the solution 
proposed was, whenever P(I,I) became negative, to 


substitute it by its positive counterpart. Accordingly 
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the following statements were introduced in the Kalman 
filter ( KLMN subroutine ) 3: 
DOe 24s Tt 1, NX 
LEP (isl). LIO.0)e PCr, =r eerie 
21 CONTINUE 
The effects of this procedure were then verified 
on the following combinations of R and 2s." 
(Gelso aes |) 
( Le], [5.J) 
Camel) nee) 
For the second case, this had as a consequence to 
hold XHAT(4) within bounds longer than before, only to 
diverge more violently towards the end ( XHAT(4)=3735. ). 
The cther two parameters had also large values. For the 
first case, the influence of the above procedure is more 
evident, since originally XHAT(4) decreased, crossed 
zero at approximately the 100th measurement ( point 
where P(4,4) became negative ) and went on decreasing. 
Vith the presence of the above statements, it behaved 
in the same way for the first 100 measurements, then 
stayed at a value of .8 until the 500th measurement, 
Searing a slow decreass through zero to reach the 
Value of -2.605 at the end . Acordingly the variance 
ft,¢), after touching zero increased to a final value 


of .242 . The same type of behaviour is also noted in 
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both XHAT(3) —— diverging with negative values before, 
with P(3,3) increasing and being constant now with 
P(3,3) increasing also, and XHAT(5). —- going away 
from zero with positive values before and going away 
from zero but with negative values and at a mech slower 
rate, with the introduction of the above statements. 
P(4,4) seemed to be bouncing between zero and some 
finite positive values, but since the velues of P were 
only recorded every 50 measurements, no valid conclusion 
can be taken about that fact. P(3,5) always increased 
and P(5,5) remained at approximately .057 . 
It was decided then to modify the structure of 

the state equations to cneck if anotner form of the 
S6quations ore converges to a solution. In Equation 
(6.15), for example the two groups forming the 
eetfictent of x(2) = u(1) could be considered as 


one parameter, being the new forms of x(4) and x(5) 


<< Gi K «= W.BG i 
ioe) = ee ene »==-8 
Se os Lx 
(7,1) 
We. BG 1 
LO X(5) mere ere genes 
ao | Ww) 
Ne, +. aes 


Anda the new form of tha state equations would then 0s 
et) =e Ot x(3) x(1) 4h x(4)) x2) =] utd) | 
eo (5) x2) ee Ole eae Aes x1) 3 


aa Oo12 ee Oley x(5)] [ x(2) ~ all) | 7 
‘> 


(ie) 


Oy x 5) eS 
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y(2) = x(1l) 
For comparison only the following valued of R and 
Po were considered : 
i) es |e 1" eval at [5] 
ii) R=[.5 ] and Po= [2s] 
Also and to check the behaviour of the system at 
high values of Q and R, the following run was made ; 
naak Q=| 1. | , R=[ 2.] ; Pe | Lal 
The main caracteristics of these runs was that the 
parameters did not oscillate, or oscillated with a very 
long period. Only in the case R-=[.2] , P=[5] , 
XHAT(3) went over 1. at the 310th measurement, 
approximately. XHAT(4) and XHAT(5) which started at 
“.o/5 and +.274 respectively, crossed each otner and 
zero at the s2Oth measurement reaching at the end the 
Values +2559 and -.422 and showing in between a 
behaviour completely linear. P(3,3) incrsased to a value 
of 5.054 at the 100th measurement and then decreassd ag 
steadily as it had increased to a final value of 6.412. 
Similerly, P(4,4) increased to a value of 5.003 in the 
Perso 15 moasurements and then decreased to a final 
value of 4.838. P(5,5) from the start decreased oscilla- 
Ting until it reached a final value of .38¢4 . This 
meowed Jike being the object of the search apart from 


the fact that the values the parameters reached at the 
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end were not to be expected on physical grounds, 
specially XHAT(5) ( this would mean that the center of 
bouyancy had to be above the center of gravity which is 
net true from the: calculated curves ). A possibility to 
@onsider is that the variables go througara cycle of 
very long period and tnerefore the convergence will take 
place later in the process. For the case R | +3] ; 
EP [ 1. | » only the value of P(5,5) decreased and the 
parameters changed their value very slowly and therefore 
nothing could be said about such a process. It could be 
an oscillatory motion of a period even longer than 
before or it could be a steady, but very slowly divergent 
processe In the last case considered, the parameters did 
not eee but the diagonal elements of the covariance 
matrix increased and quite rapidly to reach values of 
Perot , 2,259 , 48.289 , 48.195 , 52.620] . Note that 
pare tne variance of XHAT(1) and XHAT(2) was, for all 
the cases considered, simulation included, of ths order 
eieeeeo and .012 after tne first 2O to 50 measurements. 
Therefore this last case is of no consequence, 
in summary in this chapter the results from the 

identification were presented, both for the simulation 
carried out and for uss of the data. It was pointed out 
Sethe difficulty of the system to converge. A set of 


values of the different covariance matrices was selected 
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from the simulation studies as being those that offered 
more possibilities when using real data. Hovever the 
results didn't meet the expectations, being necessary 
to proceed in looking a set that would make the system 
converge - In stead it was found that some of the 
diagonal terms of P would become negative. First it was 
thought to be resultant from a mistake in the compute- 
rised scheme and later to errors introduced by the 
Semputer itself and corrected accordingly. Even then 
no solution could te found after 1000 measurements . 
Tnen another form of tne state equations was obtained 
by combining two of the parameters. This later system of 
equations presented some Interesting caracteristics, 


but nothing could be said in definite. 
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CHAPTER 8 


CONCLUSIONS AND RECOMENDATIONS 
From Chapter 7 it was quite clear that convergence 
of the estimate of the system to its true value was, in 
general, not possible to obtain. Reasons that could 
account for this failure in convergence were then pointed 
ome as follows : 
1 = Possible mistakes in the computerized 
scheme. 
2 - Not accurate prediction of the initial 
conditions<of:the' system. 
5 - Data that measured the dynamical behaviour 
was not accurate enough. 
4 - The input to the model didn't really have 
the same caracteristics ag the input to the system 
Pi tne experiments. 
To those reasons above another may be pointed 
out that could be partly responsible for the non- 
convergence of the system, namely : 
5 «- Incompatibility between the system 
structure and the particular technique that was chosen . 
With respect GO CUNO VEE St re aeOn, ho taco, 
mentioned already that this same procedure was used 
before in the identification of linear static systems 


successfully, although not a sufficient reason to 








Sefer 


completely exclude this possibility, at least decreases 
the probability of such an event. Thersfore it is recom 
mended that if future use is intended to such a form- 
lation, the program be rechecked for mistakes. 

The not accurate prediction of the inreural condi-— 
tions can be a reason of some standing if the methods 
used for prediction fall under the category of empirical 

Se semiempirical formulations from which the final 
result have not the reliability required. For the calcu- 
muon of the motion parameters for roll, for example, 
the three categories in which Chapter 3 was divided 
also provide a good basis for classification of 
Mecuracy, starting with the hydrostatic coerficients 
at the Eoeeana the non-linear coefficients at the 

bottom of the scale. It is believed that for the case 
Simeeeoi tation due to regular waves a formulation as 
meccontved by Ref. 18 and 19 for linear added mass, 
added moment of inertia and damping coefficient is 
completely adequate in the sense that the system will 
converge if the values of those parameters constitute 
mae only apemeecimation involveds  =imis cinclusiom is 
reached with information on the simulations carried 
out in Chapter 7, 

Vhen trying to digitize the data from the towing 


tank experiments, it was believed that this item 
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played a major role in the success or failure of the 
identification. This is not so anymore, when one compares 
the data from the tank with the measurements generated 

in the simulation, where the " level "‘of noise is mich 
higher than in the former case. This is not the same as 
saying that it is not necessary to care about the 
experiments. If nothing else a careful planning will 
save a lot of time and labor, waisted otherwise trying 
to correct past mistakes. 

The accurate determination of the input to the 
system as a function of time seems to play an important 
role in the success of the identification. Intuitively 
if the behaviour of two systems istto be:compared, it 

is necessary that the input be the same. Unfortunately 
this was not the case. Both methods presented in 

Chapters 3 and 4 ( using the Haskind Relations or 
assuming the major part of the excitation to be due to 
changes in bouyancy ‘dueito:the passage ofthe wave ), 
will not give very accurate results, the first because 
Lt doses not present the inataneous values of the input, 
but an expression that assumes before hand an harmonic 
tims dependency, the second because, altnougn 14 
expresses the input as the instantaneous function of 
a quantity being measured in the experiment, exprimes 


a very rough approximation of such dependence. However 





Oss 


the second was used im the Identification proccaurmom ind 
it is still believed that, por this particular case the 
Sac onGsyelds™ becren Nests. 
The incompatibility between the system structure 

and the method of identification exists because some of 
the assumptions set forth in the derivation of this 
particular method of identification were not met by the 
system. One of such assumptions was that the noise in 
the system, both environmental and measurement noise 
were uncorrelated with each other and with the state 
of the system ( Equations 6.5 and 6.6 ). This last 

Bo mt was not met by the vehicle structure since, for 
example, finite water depth effects, interference 
between the vehicle and the walls, etc., would in one 
way or another be correlated with motions of the 
vehicle. This fact as a reason for failure to achieve 
ca@nvergence gained importance when comparing the 
results of the simulation with the results of the 
identification, where the former, being bad, were a 
bit more encouraging than the last. 

With the above in mind, it seems that in order to 
improve the application of the identification techniques 
to problems of the type discussed herein, it will be 
necessary 3 


1 - Generate a vehicle environment that can be 





better approximated as an uncorrelated random process. 
This sugests the possibility of using random seas to 
generate the data for the identification. 

2 - Use another method of system identification, 
such as the model refference technique briefly 
described in Chapter 6, which do not make such restric=- 


tive assumptions. 
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HYDROSTATIC COEFFICIENTS 


1 - DISCRIFPTION OF THE COMPUTER PROGRAMS TO OB-~ 
TAIN A POLYNOMTAL SURFACE REPRESENTATION 
2 - VALUES OF THE HYDROSTATIC COEFPICIENTS FOR A 


MODEL OF THE " MARINER " HULL FORM 
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VALUES OF THE HYDROSTATIC COEFFICIENTS 
Since the identification of the motion parameters in 

coupled heave and roll was not performed, it was necessary 
to calculate the hydrostatic coupling coefficients. To cal- 
culate only the coefficients a 5 7b Kg and Kgs a much 
simpler method than resorting to the Kerwin polynomial re- 
presentation and that will probabilly yeld more accurate 
results can be used . The reason for such is that in com~ 
puterized solutions using least squares, it is not usual- 
ly possible to specify some a priori knowledge that one 
might have about the surface or the curve. For example in 


roll K(¢=0) and Kq2 are zero due to port and starboard* 


7 
symmetry. Also esl gee ~ YerGM . A least square method gave 
nonzero values fOr Une first GWor cect tl icacites. 

Looking just for the secon non zero term in the Tay- 
zor's expansion and knowing the first it's much easier to 
calculate manually the missing term. Using this procedure 
and being, for roll K, - W*GM -2.159 LBertT, the value that 
will approximate the moment curve vetter will be ; 


Between Pet.175 1/3'3K55 = 64814 LEFF? 


j) 


Pzt 350 -1.254 LBP? 


y 


p=+ 4525 -3.007 LB*FT 
Since the roll angle will go over .175 but does not 
Beech .500 in the axveriments the value choosen was: 
1/3'Ko3 = 6.814 LEXFT 
Similarly for heave the values found were 


Za = - P¥R¥ALY = ~194.709 LE/FT 


( 


2" 252 ~ 84.528 LB/FTA 
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APPENDIX 2 


LINEAR HYDRODYNAMIC COEFFICIENTS 


1 -— DISCRIPTION OF THE COMPUTER PROGRAMS TO CAL- 
CULATE ADDED MASS, ADDED MOMENT OF INERTIA AND DAMPING 
COBFPICIENTS FOR TWO DIMENSTONAL CYLINDERS IN HAVE, SWAY 
AND ROLI, | 

2 - VALUES OF ADDED MASS, ADDED MOMENT OF INERTIA 

AND DAMPING COEFFICIENTS FOR A MODEL OF THE " MARINER " 


HULL FORM IN HSAV® AND ROLL 
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PROGRAM INPUT 
Cards 1-4 FORMAT (20A4) 

The purpose of these cards is to provide titles for 
the computed output. Vhen titles are not desired blank 
cards should be provided. 

Card 5 -NOK,MAXB, NH FORMAT (316) 
NOK - Total mumber of frequencies, at which the cal- 
culations will be performed.(MAX 50) 

MAXB -Integer correspondent to the horizontal plane 
of maximum beam. Start the count at the keel. 

NH - Number of depths at which the calculations will 

be performed.(MAX 10) 
Card 6 = CAY(K),K 1,NOK HORWAT (SH 2. 7) 


CAY - Nondimensional frequencies ~~ a/g 


radian frequency 


99 
) 


half the beam for e floating cylinder 


half the draft for a submerged cylinder 
g - gravitational acceleration 
Card 7 -(DEP(J),J 1,NH),CR FORMAT(5F12.7) 

DEP - Values of the deptns of submergence, measured 
meomeuie surface to the top of the cylinder. For a floa- 
Pame cylinder, NH 1 and DEP O . 

Card 8 - MD,NON FORMAT (216) 

MD - Control integer - 1 for a submerged body 

2 for a floating body 
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NON - Number of segments in which the right hand side 
half contour is divided (MAX 45) | 
CR - Positive vertical distance from the center of 
mOotTation to the free surface, for a floating cyl i ndernwer 
heom vne center of rotation GO the inbersec Pionsoie ye 
section contour with the plane of symnetry( on the side 
closer to the surface. 
Card 9 - X(J),J 1, (NON 1) FORMAT (5F12.7) 
X(J) - Horizontal offsets of the end points of the 
line segments, measured from the vertical planr of symmes 
try to the right hand side of the contour, starting at the 
keel. The last point is the intersection of the section 
contour either with the free surface ( floating cylinder ) 
or with the vertical plane of symmetry( submerged cylinder 
Card 10 - ¥(J),J 1, (NON 1) FORMAT (5F12.7) 
Y(J) - Negative vertical offsets of the end points of 
the line segments, measured from th free surface, for a 
floating cylinder or from the intersection of the contour 
with the vertical plane of symmetry, for a submerged cy- 
linder, starting from the keel. Again only the right hand 
side of the contour will be considered. 
Card 11 - MODE FORMAT(I6) 
MODE ~ Control integer ~ 1] - heave 
2 - sway 
oO. Ss cro It 


Card 12 - FORMAT (20A4) 
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Provides a title for the ontputeol the proc are 
the title is not desired, a blank card should be provided. 
If more than one mode is desired, adicional pairs of 
@ords 1] and le should be provided, accordingly. 
Card 13 - MODE FORMAT (I6) 
MODE - Control variable 
O - End of data 
Negative integer = A new geometry follows 


Starting Card 5. 
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APPENDIX Ss 


NON-LINEAR MOTION PARAMETERS 


1 - STERMINATION OF HYDRODYNAMIC COEFFICIENTS 


USING EMPIRICAL METHODS 


’ 
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L - ADDED MASS AND ADDED MOMENT OF INERTIA 


Using Ref. 20 pemolmavowoe - Heave 


7 ae 2 
am = 4% (0% "2 ] Zz y- (dx) (A3.1) 
whe re - water density 1.94 LBxsEc’/rr4 
~ lenght of the model 5.5 FT 


coefficient from ite. 7 se. 


ic 
$ 
My 


- contour of the waterline, aproximated by a 


parabola of tne form 


n 
y- B/2| 1 ~ +5 (A362) 
B -» beam of the model Pais Nk 
fiers Cay (1-0eN) 
Cw - waterplane coefficient - Ay / BEL = 724 


ra = coefficient, calculated acceovecane to 


4, > 1/ | 1+.8*B/L | = +395 (A3.3) 


the following form may be obtained 
Am = o772xFx @ 3 “LeB*: ne / (Le C, ae oes) 


Pimeoi ng in the values 
Am=1.41 LBYSEC"/FT 
eS 
Being the mass of the vehicle += m 1.43 LB*SEC /FT 


this result may be checked by measuring the natural period 








== 


of heaving motions and using results from the hydrostatic 


calculations Z,=- *g*A , = -194.05 LB/FT 
oh 
2 
- 2, 
w= = 98.20 see 2) 
m+Aam 


The value found for the natural undamped frequency was 
8.038 (SEC-1). 
Similarly for roll, the added moment of inertia is 
expressed as a percentage of the moment of inertia. This 
is found, from data reffering to full scale vehicles, acor-~ 


ding to the formula 


c 


2 2 
i. = w/a ( (BY xcw)/11.48C, + pD~/12 ) TON“SEC (A3.5) 


i ©02855 LBSEC” 
From fige 24 of the above refference 


a Iyx ae 28 OI» 


There fore 


.03665 LBYSEC“ 


M 


I Pe eull 
KX 2, 89.4 


From the hydrostatic calculations 
Kg=- WatGH = -2.159 LBSFT 
The resulting natural frequency is then 


GO, =. 7.7 Sue(-1) 


while the value found in the measurements was 


Won = 4.14 ne (el) 








The values of added mass and of addes moment of 
inertia, calculated in Appendix 2 were : 
A(3,3) = 1.24 LB«sEc®/rr 
A(4,4) = .0176 LB*SEC“sFr 
The value of roll added moment of inertia is seen to 


be quite different the results of both methods are compared. 


z~ = DAMPING COEFFICIENT 
For heave the following expression, taken from page 


153 of the above refference was used : 
W 


Z 1/20 2YAMD-B.G,) © (B/A) « (E6E/0 (ASS) 


where the parameters are as before and also ;: 


» « wave lenght = 3.11 FT for w=8.04 SEC(-1) 


Xx" 2 coefficient taken from (20) fig. 52 = 4460 
T « aes Of the vehicle = 322 71 

:&, = coefficient taken from (20)fig. 53 = 200 
For the values above : Z_= -2.336 LE“S&C/FT 


while when calculated by the methods of Appendix 2 

Quy = 704543 LBXSEC/FT which is three times the 
Value presented above. 

For roll, the calculation of linear and non-linear 
damping coefficients were carried simultaneously, for each 

mreansverse section of the vehicle, using Equations Zl, 22, 

2o or 24 of page 148 of the above refference, deppending on 
the shape of the section. The value of the damping for the 


whole vehicle was obtained by summing the sectional results 
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over the lenght of the vehicle ina" strip theory " waye 


Therefore, using the following equations the damping is 


calculated. 
2 
Kp, ~1/2 eb, a1, Con (5, B/E) [PY es ee SCT) 
378 


72(B,/T, bia 3.07)+11 (A368) 
p 


1,:4 ; ae ee ae 
Kp, ~2fP¢414 |CQ3 (7Q5B,/T,) [phen ae 


Di 





es cea oe oe 
Kos wep, ol, Pa ee By Via Lary oe ~ geo 5 
2 
2b ‘a1 57 |p , Bibs/T72-82) 4 (AS .10) 
Sos (Pat | OTP FE TO 





See INET 


where Eq. A3.7 and A3.8 are used for rectangular sections 
of different value of ro ry/b; , Eqs AS.9 is used with 
elliptical sections and Eq AS.10 with terminal section of 


shape close to triangular, the following values were 


~ 
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found for linear and non linear damping coefficients: 


Ky = -.00786 LBFTSEC ~ W= 4.145.SEC(-1)} 
— -,.06835 LB“FTsSEC - w= 8.290.SEC(-1) 
KD ip - O1LEO5 LBS FT+SEC“ - all frequencies 


The values fourd in Apenndix 3 for the linear damping 
coefficients were respectively -.0022 LB“ FT“<SEC 
-.0363 LB*FTsSEC 
For the nonlinear heave damping coefficient, the va-= 
lues for each section were obtained by considering each se- 
ction substituted by an equivalent cylinder of the same 
diameter when the model moves away from the water ( assu- 
ming the only drag present to be separation drag ) and sub- 
stituted by a cylinder of the same lenght of wetted con-~ 
tour when the model moves into the water ( assuming the on- 
ly drag present to be skin friction drag ) . Data fron 
infinite cylinders was then used, the correspondent values 
of drag coefficient averaged over a cycle, and the results 
summed over the lenght of the model. The velocity entered 
in the calculation of the Reynolds Number was a RNS type 
of velocity taken from the records of the experiments. 
The final value obtained was: 


Re aie 
Zu |yl = —~L-767 LEsSEC™/FI 
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APPENDIX 4 


CIRCUITRY AND DIGITIZING PROGRAM 


Ll - PROGRAM TO DIGITIZE DATA AND CORRESPONDENT 
PATCHING IN THE ANALOG COMPUTER 
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APPENDIX 5 


IDENTIFICATION PROCEDURE 


L - DISCRIPTION OF THE COMPUTER PPOGRAMS USED IN 
THE IDENTIFICATION PROCEDURE 
2 - EXAMPLE OF THE BEHAVIOUR OF VARIABLES AND RE- 


LATED VARTANCES DURING IDENTIFICATION 
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